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Abstract. We investigate distortions in the velocity fields of disc galaxies and their use to reveal the dynam- 
ical state of interacting galaxies at different redshift. For that purpose, we model disc galaxies in combined 
N-body/hydrodynamic simulations. 2D velocity fields of the gas are extracted from these simulations which we 
place at different redshifts from z—Q to z=l to investigate resolution effects on the properties of the velocity field. 
To quantify the structure of the velocity field we also perform a kinemetry analysis. If the galaxy is undisturbed 
we find that the rotation curve extracted from the 2D field agrees well with long-slit rotation curves. This is not 
true for interacting systems, as the kinematic axis is not well defined and does in general not coincide with the 
photometric axis of the system. For large (Milky way type) galaxies we find that distortions are still visible at 
intermediate redshifts but partly smeared out. Thus a careful analysis of the velocity field is necessary before using 
it for a TuUy-Fisher study. For small galaxies (disc scale length 2 kpc) even strong distortions are not visible in 
the velocity field at z ~ 0.5 with currently available angular resolution. Therefore we conclude that current distant 
TuUy-Fisher studies cannot give reliable results for low-mass systems. Additionally to these studies we confirm 
the power of near-infrared integral field spectrometers in combination with adaptive optics (such as SINFONI) to 
study velocity fields of galaxies at high redshift (z~2). 

Key words. Galaxies: kinematics and dynamics - Galaxies: interactions - Methods: numerical 



1. Introduction 

Recently it has become technically feasible to observe the 
full 2D velocity field (VF) of local galaxies in optical wave- 
bands using integral field units (IFUs) such as SAURON 
(e.g. Ganda et al. 2006) or Fabry-Perot interferometry 
(e.g. Chemin et al. 2006, Garrido et al. [2002)) . For inter- 
mediate and high redshift galaxies, however, there are by 
now hardly any observational studies of 2D velocity fields 
available. 

Flores et al. (2006) observed the 2D velocity field of 
35 galaxies at intermediate redshift (0.4 < z < 0.75) using 
FLAMES/GIRAFFE at VLT. One aspect was to investi- 
gate the redshift evolution of the TuUy-Fisher relation. A 
different approach, used by our group, was presented by 
Ziegler et al. (2006) and Kutdemir et al. (2007) who uti- 
lize multiple-object spectroscopy from the VLT with dif- 
ferent slit positions on each galaxy in order to construct 
the full velocity field for each galaxy. Most of the other 
studies of distant, faint, and small galaxies are still based 
on slit spectroscopy (e.g. Vogt 2001, Bohm et al. 2004). 
To account for distortions and irregularities in the veloc- 
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ity fields is in both cases critical, especially when aiming 
at a distant Tully-Fisher study. In two recent papers we 
showed that observational constraints and galaxy-galaxy 
interactions can severely influence the determination of 
the rotation curve of observed disc galaxies (Kapferer et 
al. 2006, Kronberger et al. 2006). 

In this paper we investigate to what extent the full 
2D velocity field of a galaxy can be used to gain infor- 
mation on its internal kinematics and can hence improve 
the quality of e.g. Tully-Fisher studies. This question is 
also important to possibly disentangle different interac- 
tion processes by mapping 2D velocity fields and to study 
their impact on galaxy evolution. We focus on the question 
how the visibility of distortions depends on the redshift of 
the observed galaxy, i.e. the actual spatial resolution of the 
galaxy. For that investigation we place a galaxy at differ- 
ent redshifts and bin the velocity according to the spatial 
resolution at this redshift. A similar study for observed 
galaxies was presented by Epinat et al. (2006). 

Additionally to the examination of 2D velocity fields 
of intermediate redshift galaxies from optical spectroscopy 
we investigate the performance of near-infrared integral 
field spectrographs that are used together with adap- 
tive optics. As a prototype we take the characteristics of 
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SINFONI at the Very Large Telescope of the European 
Southern Observatory to study the velocity fields of galax- 
ies at z~2. For example, Genzel et al. (2006) observed the 
velocity field of a massive protodisc at z=2.38 detecting 
an ordered rotation without any hint for a major merger 
event in the system. 

Very recently also Jesseit et al. (2007) analysed 2D 
velocity fields of simulated galaxies. They focused on sim- 
ulated disc merger remnants and found that many differ- 
ent kinematical phenomena can be observed in the stel- 
lar velocity maps, such as kinematic misaligned discs or 
counter-rotating-cores. For the analysis they also used the 
kinemetric method of Krajnovic et al. (2006) as we do in 
this paper. They did, however, not investigate the redshift 
dependence of the 2D velocity fields. 

The paper is organised as follows. In Sect. 2 and 3, 
we describe the simulations, the interaction geometries, 
and the way we extract realistic 2D velocity fields from 
the numerical data. In Sect. 4 the results for different in- 
teraction scenarios and their dependence on the angular 
resolution are presented. We end with a summary of the 
main conclusions in Sect. 5. 

2. Simulations 

In this work we use some of the simulated systems pre- 
sented in Kapferer et al. (2005), which were subsequently 
used by Kapferer et al. (2006) and Kronberger et al. 
(2006). The simulations were carried out with the N- 
body/SPH code GADGET-2 developed by V. Springel 
(see Springel 2005 for details). In this code the gas of 
the galaxies is treated hydrodynamically and prescriptions 
for cooling, star formation, stellar feedback, and galac- 
tic winds are included (Springel & Hernquist, 2003). The 
coUisionless dynamics of the dark matter and the stellar 
component is calculated using an N-body technique. The 
initial conditions were built according to Springel et al. 
(2005), based on the analytic work of Mo et al. (1998). 
The two model galaxies were chosen such that they repre- 
sent a Milky Way type and a small spiral galaxy, with the 
mass ratio of the two galaxies being 8:1. Therefore, the 
total mass of the model galaxies A and B is 1.34 x 10^^ 
h-^ Mq and 1.67188x10" h''^ Mq, respectively. The 
combined N-body/SPH simulation then calculates 5 Gyr 
of evolution. For every time step, we know the velocity of 
each particle and can hence extract realistic 2D velocity 
fields of the gas. 

Concerning the spatial alignment and impact parame- 
ters, we follow the notation introduced by Due et al. (2000) 
to describe the interaction geometry. The parameter b 
corresponds to the minimum separation of the galaxies' 
trajectories, as if they were point masses on Keplerian or- 
bits. Additionally, two angles, 8 and $ define the spatial 
orientation of the disc. 

For the complete sample, we selected the alignments 
in such a way as to cover as many geometries as possible, 
including minor and major mergers and fly-bys (achieved 
by increasing the minimum separation). For this analysis 



Table 3. Physical resolution according to 1 arcsec angular 
resolution for different redshifts in the standard ACDM 
cosmology. 



Redsliift physical resolution according to 1" [kpc] 



0.05 


0.98 


0.1 


1.85 


0.2 


3.30 


0.3 


4.46 


0.4 


5.38 


0.5 


6.11 


0.8 


7.51 


1.0 


8.01 


2.0 


8.38 



we just use three different interaction geometries, which 
we list in Table [TJ These three simulations allow us to 
study different classes of kinematical distortions, detailed 
in Sect. [3l 

Additional important quantities of the simulations, as 
e.g. the particle numbers, are summarized in Table [H 
Throughout the paper we adopt the standard ACDM cos- 
mology with r^A = 0.7, = 0.3, and h=0.7. 

3. Extraction of realistic 2D velocity fields 

In order to construct realistic 2D velocity fields, we project 
all gas particles of the N-body/SPH simulations onto a 
Cartesian, equally spaced grid. The spacing is chosen such, 
that the spatial resolution at the assumed redshift of the 
model galaxy corresponds to the angular resolution of 
current state-of-the-art observations. For the whole inves- 
tigated redshift range from z— to z=1.0 we adopt an 
angular resolution typical for IFU or FPI observations, 
namely 0.4" (e.g. Chemin et al. 2006). The angular res- 
olution of SAURON, for example, would be 0.3" or 0.9". 
For intermediate redshifts we additionally chose 0.25", 
which is the pixel scale along the F0RS2 slit as used by 
Ziegler et al. (2006) and 0.52", which is the lenslet size of 
FLAMES/GIRAFFE at VLT (e.g. Flores et al. 2006). We 
calculate for each redshift the physical resolution accord- 
ing to the given angular resolution using the concordance 
cosmological model (see Table [3]). The velocity field of the 
galaxy is binned using this spatial resolution. 

From the knowledge of the full 3D velocity fields and 
the interaction history of the galaxies in our simulations 
we define three kinematical classes, in principle analogue 
to classifications in Flores et al. (2006) or Krajnovic et al. 
(2006): 

— Undisturbed Rotation: The galaxy is not interact- 
ing and the rotation curve has the "classical" shape, 
rising in the inner part and turning over to a flat 
regime. 

— Disturbed Rotation: Rotation is still predominant 
but disturbed by a minor interaction, e.g. a minor 
merger event. Non-rotational components of the ve- 
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Table 1. Interaction parameters for the simulations. C2: the two interacting model galaxies. Model galaxy A is a 
Milky Way type spiral galaxy while model galaxy B represents a small spiral galaxy, with the mass ratio of the two 
galaxies being 8:1.; C7: minimum separation [kpc]; C8: initial relative velocities of the two galaxies [km/s].In column 
C9 we give the initial orbital energy of the encounter in Joule. 



Simulation 
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62 


C7 


C8 
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50 


250 


4.8 X 10"" 
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120 


9.7 X 10^° 
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4.8 X 10^^ 
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A-B 
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250 


4.8 X 10^^ 
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A-A 
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350 


1.1 X 10^^ 



Table 2. Particle numbers, mass resolution and gravitational softening used for the two model galaxies. Additionally 
the circular velocity of the halo at r2oo and the disc scale length of the initial conditions are given. 





particle number 


mass resolution 


softening length 


halo circular velocity 


radial disc 






[h~^ Mq /particle] 


[h"^ kpc] 


[km/s] 


scale length [hr'^ kpc] 


Galaxy A: 


Dark matter halo 


30000 


4.2 X 10" 


0.4 


160 




Disk coUisionless 


20000 


3.3 X 10^ 


0.1 






Gas in disk 


35000 


8.4 X lO"* 


0.1 




4.5 


Galaxy B: 


Dark matter halo 


30000 


5.3 X 10'' 


0.4 


80 




Disk coUisionless 


20000 


4.2 X 10'* 


0.1 






Gas in disk 


35000 


1.1 X lO"* 


0.1 




2.25 



locity field are also clearly visible as peculiarities in 
ID rotation curves (Kronberger et al. 2006). 
— Distorted Velocity Field: The velocity field is heav- 
ily distorted by a strong merger event, i.e. an interac- 
tion with a galaxy of similar or higher mass. It shows 
a very complex pattern and has no regular rotation 
curve anymore. 

The main question that we investigate in this paper 
is how such classifications depend on the redshift of the 
observed galaxy, i.e. the actual spatial resolution of the 
galaxy. In this context the seeing plays a crucial role, as 
it typically exceeds the angular resolution of the instru- 
ment. To simulate seeing effects on our velocity-field mea- 
surements, a convolution with a Gaussian point spread 
function was applied. We adopted a value of 0.8" for the 
FWHM of the Gaussian seeing, which is a typical value for 
ground based observations (see e.g. Jager et al. 2004). The 
appearance of the velocity fields is therefore dominated by 
the seeing. Note that we do not calculate an evolution of 
properties of the galaxy with redshift but study how a 
given kinematical state of a galaxy is observed at different 
redshifts. 

All studies mentioned above are mainly seeing dom- 
inated. We additionally study the possibility to identify 
mergers at high redshift (z'^2) with an adaptive optics in- 
strument such as SINFONI at the VLT. For this investi- 
gation we adopt an angular resolution of 0.15" as achieved 
by Genzel et al. (2006). 



4. Results 

In the subsequent sub-sections we will systematically in- 
vestigate velocity fields from each of the kinematical 
classes defined in Sect.|3l The focus lies on the visibility of 
kinematical distortions as a function of redshift. In order 
to quantify the distortions and to interpret the partly com- 
plex structures in the velocity fields we use the kinemetry 
package of Krajnovic et al. (2006). This analysis is based 
on the assumption, that the mean velocity along best fit- 
ting ellipses can be reproduced by a cosine law, i.e. 

V{a,-^)^VQ + Vc{a)cos^, (1) 

where a is the length of semi-major axis of the ellipse, ^' 
is for discs the azimuthal angle measured from the ma- 
jor axis in the plane of the galaxy. Note that in general 
a simple sine correction for the inclination of the disc is 
applied in observations too. The velocity fields presented 
here have not been corrected in this way as all of them 
were extracted at the same inclination i=35° (except for 
those presented in Sect. 14. 4p . 

The position angle F and the axial ratio {q — h/a) 
of the ellipses are calculated as a function of radius from 
the galactic centre. Deviations from that cosine law are 
measured using an harmonic expansion along the ellipses, 
i.e. 

N 

V{a, ^) = Y, kn{a)cos[n{^ - 0„(a))], (2) 
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where (f) is the phase coefficient. Higher order Fourier 
terms and radial changes of F or g are used to quantize de- 
viations in the velocity field from a simple rotation. There 
is an intrinsic relation of the surface brightness of the 
galaxy and the line-of-sight velocity. Both are moments 
of the underlying distribution function, where the surface 
brightness is the zeroth's order moment and the line-of- 
sight velocity is the first order moment. Note, however, 
that for gas, as in our case, the situation is different from 
collisionless stellar systems, where the distribution func- 
tion for simplified systems can be calculated analytically. 
As we deal here mostly with rotating disc-like structures 
the use of ellipses is a natural choice. For distorted and 
elliptical systems this assumption will be violated and the 
deviation is quantified by the Fourier analysis. In prin- 
ciple the calculation of the best fitting ellipses works in 
a two step approach as follows. First, higher order terms 
are minimized, which do, according to their analysis, not 
carry the major contribution to the map. This happens 
on a grid of position angles (F) and flattenings (q) and 
leads to initial values (Tmin) and {qmin) for a second fit. 
In this second fit all parameters of the ellipse are taken 
into account and the Fourier analysis can be performed 
to arbitrary order. For more details on the algorithm see 
Krajnovic et al. (2006). 

4.1. Undisturbed velocity fields 

As a first step we investigate the simple case of an undis- 
turbed spiral galaxy, i.e. an example for the first kine- 
matical class presented in Sect. [S] This class is partic- 
ularly important for TuUy - Fisher studies, as only for 
these galaxies, the maximum rotational velocity Ymax is 
an approximation for the virial mass of the system. For 
this investigation we use galaxy A in simulation 1 (see 
Table [T]), at a snapshot when the two galaxies are still 
well separated. The adopted inclination of the galaxy is 
35° (90° is defined to be edge-on). The velocity field is 
very regular at all redshifts. The appearance of our VFs 
is dominated by seeing, as we adopted a constant value of 
0.8" for the FWHM, which is always larger than the an- 
gular resolution we use. Both effects, the worse sampling 
at higher redshifts and the large seeing smear out the ve- 
locity field when shifted to higher z. In Fig. [T] we present 
the VFs of an undisturbed disc galaxy as seen at two dif- 
ferent redshifts (z=0.1 and z=0.5). Overlayed on the VF 
are the best fitting ellipses from the kinemetry analysis. 
At a redshift of 0.1 there is some structure visible in the 
VF, e.g. a small twist of the position angle F towards the 
centre. Also in the disc some fluctuations of the rotational 
velocity are present. These small structures in the VF are 
completely smeared out at redshift z — 0.5. 

Quantitatively the differences can be seen in Fig. [2] in 
the radial profiles of the kinemetric properties, calculated 
using the kinemetry programme. The program also calcu- 
lates formal la errors from the covariance matrix, which 
are shown as error bars in the plots. These uncertainties 




Fig. 1. 2D velocity field of an undisturbed modelled disc 
galaxy as seen at redshift 0.1 (left) and at 0.5 (right). 
Overlayed as contours are the best fitting ellipses from 
the kinemetry analysis. 



are estimated from the measurement uncertainties in the 
kinematic data. As we do not account for the spectral reso- 
lution of the instruments we use the scatter of the velocity 
field. Note that this is just a formal uncertainty, which in 
reality is higher, also because systematical errors add up. 

While at z = 0.1 both, F and q show variations with 
radius, they are almost constant at z — 0.5. The first or- 
der moment fci corresponds to the rotation curve (RC) 
of a spiral galaxy or more generally spoken to the bulk 
motion in the velocity field. For the undisturbed velocity 
field presented here it shows the typical behaviour of a 
rotation curve, i.e. rising in the inner part and turning 
over to a fiat regime. This undisturbed shape is present 
at redshift z = 0.1 and at z — 0.5. In Fig. [3] we show the 
rotation curves of the system at redshift z = 0.5 obtained 
from the 2D velocity field (triangles) and from a simulated 
slit (asterisks), extracted as described in Kronberger et al. 
(2006). The modelled slit was placed over the major axis 
of the galaxy with a slit width of 1" . The RCs extracted in 
these different ways agree very well, as expected for such a 
regular velocity field. As a dashed line we plot in the same 
figure the RC extracted at redshift z = 0.1, which agrees 
reasonably well with the higher redshift RCs. This demon- 
strates the principle power of 2D velocity fields for distant 
TuUy-Fisher studies if regular, undisturbed velocity fields 
are considered. 

The last row in Fig.[2]shows a quantitative measure for 
distortions in the 2D velocity field. The fifth order term of 
the harmonic expansion represents complex, kinemati- 
cally separate components in the velocity field. In the plot 
the ratio k^/ki is shown. For the regular VF presented in 
this section, this ratio is small, generally below 0.1 for 
all redshifts. However, this value is very sensitive to the 
resolution of the velocity field. While at z = 0.1 there 
are some radii with increased k^/ki, these signatures are 
completely smeared out at higher redshifts. At z = 0.5 
the ratio is almost constant below a value of 0.02. This 
already indicates problems with identifying distortions in 
2D velocity fields at higher redshifts. This issue will be 
addressed in the next section. 
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Fig. 2. Radial profiles of the kinemetric properties, calcu- 
lated using the kinemetry programme for an undisturbed 
Milky Way type galaxy at redshift z — 0.1 (top) and 
2 = 0.5 (bottom). The position angle T and the flattening 
q of the best fitting ellipses as well as the first and the fifth 
order Fourier terms fci and are plotted as a function of 
radius. 



At redshift z = 0.5 we study the effects of angular 
resolution on VF properties by considering in addition to 
the standard angular resolution of 0.4" two more values. 
We adopt additionally the 0.25" of Ziegler et al. (2006) 
in their sample using F0RS2 spectroscopy and 0.52" as 
for FLAMES/GIRAFFE at VLT (e.g. Flores et al. 2006). 
At the FLAMES/GIRAFFE angular resolution the VF is 
even more smeared out than at 0.4" resolution. However, 
the radial behaviour of the kinemetry parameters is only 
slightly affected. Especially the RCs and their maximum 
value Yrnax agree reasonably well for all three angular res- 
olutions. Again, with respect to subsequent sections, we 
point out that this holds only for regular velocity fields. 

4.2. Disturbed rotation 

In Sect. [3] we have defined the kinematical class 'disturbed 
rotation' as a rotating system with a clearly visible distor- 
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Fig. 3. Rotation curves of the isolated galaxy extracted 
from the 2D velocity field (triangles) and using a modelled 
slit (asterisks) at redshift z=0.5. The dashed line shows 
the corresponding RC at redshift z=0.1 obtained from a 
modelled slit. 



tion. As an example for this class we investigate a model 
galaxy undergoing a minor merger (simulation 2). At this 
snapshot the distance between the two galaxies is smaller 
than the radius of galaxy A, i.e. galaxy B encounters the 
Milky Way type galaxy A in its first passage. The small 
galaxy permeating the more massive galaxy A is clearly 
visible in the velocity field, as irregular structure in the 
lower right part of the VF presented in Fig. |31 At red- 
shift z = 0.1 the kinemetry analysis shows strong vari- 
ations of the position angle F and the flattening q with 
radius (see Fig. [5]). The radial fluctuations are stronger 
than for the isolated galaxy. The first order Fourier term 
fci, which corresponds to the rotation curve also shows 
distortions and the ratio k^/ki is also significantly higher 
than in the isolated case. Note that the shape and the 
maximum value of the RC are completely different from 
the regular RC. Therefore, this galaxy could not be taken 
for TuUy-Fisher studies. We note, that also RCs extracted 
from this galaxy with a slit show clear signatures of a dis- 
tortion (Kronberger et al., 2006), especially an asymmetric 
shape with a rising side pointing towards the interaction. 
However, with the help of 2D velocity fields the nature 
of the interaction becomes a lot more accessible. The al- 
most undistorted left hand side of the VF points towards 
a very asymmetric interaction, which suggests a tidal in- 
teraction. This asymmetry is less visible at z = 0.5, where 
the distortion is smeared out. Simple eye-balling might 
even lead to a misclassification as 'undisturbed rotation'. 
However, by a more detailed comparison with Fig. [T] some 
differences are clearly visible. The most striking difference 
is the presence of a clear kinematic axis in the regular VF, 
which is not present in Fig. 3) 

We present again the rotation curve extracted from the 
2D velocity field and using a simulated slit in Fig. [6l In 



6 



T. Kronberger et al.: 2D velocity fields of simulated interacting disc galaxies 




Fig. 4. 2D velocity field of a modelled Milky Way type disc 
galaxy disturbed by an ongoing minor merger event (i.e. 
the small galaxy B is permeating galaxy A in its first pas- 
sage) as seen at redshift 0.1 (top left) and 0.5 (top right). 
Overlayed as contours are the best fitting ellipses from 
the kinemetry analysis. In the lower panel the projected 
gas distribution at this snapshot is shown to illustrate the 
interaction geometry. 



contrast to the case of the isolated galaxy, the RCs differ 
significantly in the interacting system. The reason for this 
difference is the fact that in an interacting system the 
kinematic axis is not well defined and does in general not 
coincide with the photometric axis of the system. Thus, 
the slit does in general not follow the kinematic axis of 
the galaxy. Only the knowledge of the full 2D velocity 
field allows to determine the kinematic axis. 

The plots of the radial behaviour of the kineme- 
try parameters in Fig. [5] do not yield clear evidence 
for a distorted VF at z=0.5. Only the outermost data 
point of fcs/Zci shows a significant increase. The situ- 
ation is even worse for the resolution obtained with 
FLAMES/GIRAFFE. Although the maximum rotational 
velocity of the RC becomes similar to the one of the iso- 
lated galaxy due to the smearing out of the VF, a mis- 
classification as undistorted galaxy would result in sys- 
tematic errors in TuUy-Fisher studies. The luminosity of 
the system is possibly enhanced due to a merger induced 
starburst. 

In Fig. [7| we show the 2D velocity field of this example 
for redshift z = 0.05, z = 0.3, z = 0.8, and z = 1.0. It 
is clearly visible, how substructures are smeared out at 
higher redshift. 

We studied the effects of different interaction geome- 
tries on the results presented above. For this purpose we 
considered a counter-rotating unequal mass merger (sim- 
ulation 4) and an unequal mass merger, where the disc of 





J*" nJ 

1.0 IS ajj 3.0 

Fig. 5. Radial profiles of the kinemetric properties, calcu- 
lated using the kinemetry programme for the interacting 
galaxy A of Fig. [4] (minor merger) at redshift z = 0.1 (top) 
and z = 0.5 (bottom). The position angle F and the flat- 
tening q of the best fitting ellipses as well as the first and 
the fifth order Fourier terms k\ and are plotted as a 
function of radius. 



galaxy B is initially perpendicular to the disc of galaxy 
A (simulation 5). The appearance of the velocity fields 
does of course differ for the different interaction scenar- 
ios. However, the main conclusions drawn so far do not 
change. The merger is clearly visible in the asymmetric 
shape of the velocity field, the distortions are, however, 
smeared out at higher redshift. 

Of course it is more probable to observe an interacting 
system in a pre- or post merger phase than in the rela- 
tively short period where the two galaxies permeate each 
other. Therefore, we study the same system 100 Myr after 
the direct encounter described above, i.e. after the first 
pericentre passage. This is the period between the first 
passage and the complete merging of the two galaxies. The 
small galaxy went through the gaseous disc of the massive 
galaxy and left a disturbed velocity field. In the direct 
image the interaction is not visible anymore at redshift 
z = 0.5. In Fig. [5] we show the distribution of the stellar 
mass in the galaxy as 2D image and as a radial profile. By 
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Fig. 6. Rotation curves of the interacting galaxy extracted 
from the 2D velocity field (triangles) and using a modelled 
slit (stars) at redshift z—Q.5. The solid line shows the cor- 
responding RC from the 2D field at redshift z—0.1. 




Fig. 7. Evolution of the 2D velocity field appearance with 
redshift: a) z = 0.05, h) z = 0.3 c) z = 0.8 d) z = 1.0. 
Overlayed as contours are the best fitting ellipses from the 
kinemetry analysis, z = 0.1 and z = 0.5 are shown in Fig. 

m 



assuming a certain constant mass-to-light ratio this could 
be translated to a light distribution. The profile can be 
fitted by an exponential law. The velocity field shows still 
some signatures of the interaction (see Fig [9]) at lower red- 
shift but at z=0.5 these are heavily smeared out. Thus, 
the velocity information is an important complement to 
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Fig. 8. The distribution of the stellar mass in the galaxy 
as 2D image (left) and as a radial profile (right), which 
translates under an assumption of a constant mass-to- light 
ratio to a surface brightness profile. The galaxy is shown 
100 Myr after the first pericentre passage of galaxy B at a 
resolution corresponding to an artificial redshift of 2 = 0.5. 
The exponential fit to the profile is shown as a red line. 
The physical field of view of the 2D image is 36 x 36 kpc. 
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Fig. 9. 2D velocity field of a modelled disc galaxy 100 
Myr after a minor merger event as seen at redshift 0.1 
(left) and 0.5 (right). Overlayed as contours are the best 
fitting ellipses from the kinemetry analysis. 



the morphological analysis when studying interactions of 
galaxies but only if the VF is sufficiently sampled. 

Kinematic disturbances from interactions are expected 
to fade within a few rotation cycles (< 1 Gyr) (e.g. Dale 
et al. 2001, Kronberger et al. 2006). Only 200 Myr after 
the direct encounter presented above neither the direct 
image nor the VF show clear signatures of the interaction 
anymore. In a more detailed analysis, however, signs of 
the interaction can still be found. About 1 Gyr after the 
first passage, when both galaxies are well separated but 
approaching each other again, galaxy A would be classified 
as undisturbed, even when observed at low redshift. 

The analysis so far was made for the rather large 
galaxy A, which has a radial scale length of about 4.5 kpc. 
Therefore, the velocity field is sampled by many pixels. 
The situation is completely different for smaller galaxies, 
as e.g. model galaxy B. We investigate a 'disturbed rota- 
tion' of galaxy B using a snapshot of simulation 3, an equal 
mass merger of two small galaxies B. Equal mass mergers 
massively disturb the velocity field of the gas of the in- 
teracting galaxies. Most of the gas is converted into stars 
by a merger induced starburst or lost to the intergalac- 
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Fig. 10. 2D velocity field of a model galaxy B disturbed 
by an equal mass merger, 300 Myr after the first pericentre 
passage of the two galaxies, as seen at redshift 0.1 (left) 
and 0.5 (right). 



tic medium by tidal forces (e.g. Kapferer et al., 2005). We 
choose a snapshot about 300 Myr after the first encounter, 
when again some regular rotation has established. 

In Fig. (TU] the 2D velocity field of the galaxy for this 
snapshot is shown at redshift z = 0.1 and at z = 0.5 with 
an angular resolution of 0.4". While at z = 0.1 the dis- 
tortions in the VF are clearly visible, they are completely 
smeared out at z = 0.5 leaving a regular VF, that would 
be classified as 'undisturbed rotation'. This classification 
is supported by the radial behaviour of the kinemetric 
quantities shown in Fig. [TTl although at z — 0.5 only two 
ellipses were fitted due to the small number of available 
pixels. Such undetected distortions introduce an enormous 
source of systematic errors to distant TuUy-Fisher studies. 
Therefore we conclude that current distant Tully-Fisher 
studies cannot give reliable results for low-mass systems 
if the velocity field is not sampled sufficiently. 

In Fig. [12] we present the 2D velocity field of galaxy 
B for redshift z = 0.05, z = 0.2, z = 0.3, and z = 0.4. 
Gradually the distortions are smeared out, leaving an al- 
most regular VF at z = 0.4. Also for this interaction sce- 
nario we investigate how long the 2D velocity field main- 
tains the disturbed features. After the first pericentre pas- 
sage of the two galaxies the velocity field settles again to 
an undisturbed state. As for the unequal mass merger dis- 
cussed above, the strongest features disappear after sev- 
eral hundred Myr while an undisturbed velocity field is 
again present after about 1 Gyr. Note that in this pe- 
riod the two galaxies are approaching each other again for 
their second passage, but still without any direct interac- 
tion yet. 

We also investigate a major merger between two large 
galaxies A where the velocity field spreads over more pix- 
els. In this case the situation is, as expected, slightly better 
than for the small galaxy B. Especially for the kinematical 
analysis more ellipses can be fitted even at redshift z=0.5. 
Nevertheless a similar trend is observable as for galaxy B. 
While the distortions are clearly visible at low redshift, 
they are subsequently smeared out when the galaxy is 
placed at higher redshift. In Fig.[T3]we show the 2D veloc- 
ity field and the radial profiles of the kinemetric properties 
as seen at redshift z=0.5 for model galaxy A disturbed by 





OS 10 



1.* ?fl J,s 













1 


lid 

2ca 












1 

114 












da 






















i 


U 2. 










































I 




4 


OA 











Fig. 11. Radial profiles of the kinemetric properties, cal- 
culated using the kinemetry programme for galaxy B 300 
Myr after an equal mass merger for redshift z = 0.1 (top) 
and z — 0.5 (bottom). The position angle F and the fiat- 
tening q of the best fitting ellipses as well as the first and 
the fifth order Fourier terms ki and are plotted as a 
function of radius. 



an equal mass merger, 100 Myr after the first passage of 
the two discs. The appearance of the velocity field and 
the radial behaviour of mainly the flattening q point to- 
wards an interaction (compare to Figs. [T] and [5] for the 
undisturbed case). Small scale distortions are, however, 
smeared out. 

Complementary to these investigations of VFs at in- 
termediate redshift, we study the appearance of the 2D 
velocity field at z=2 as observed with SINFONI at VLT. 
Due to the use of adaptive optics, these observations are 
not seeing limited and have a high angular resolution (we 
adopt 0.15" for our study). As a test case we take galaxy 
A during an ongoing minor merger in simulation 2. The 
velocity field of this snapshot for seeing limited observa- 
tions and various redshifts was shown in Fig.[71 The galaxy 
has a radial exponential disc scale length of 4.5 kpc and is 
therefore comparable in size to the one observed by Genzel 
et al. (2006). The velocity field of this snapshot as seen at 
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Fig. 12. Evolution of the 2D velocity field appearance 
with redshift of a model galaxy B disturbed by an equal 
mass merger, 300 Myr after the encounter: a) z = 0.05, b) 
z = 0.2 c) z = 0.3 d) z = 0.4. Additionally the same 2D 
velocity fields as seen at z = 0.1 and z — 0.5 arc shovifn in 

Fig. [in] 



redshift z=2 with an instrument such as SINFONI is pre- 
sented in Fig. [T31 Peculiarities in the VF caused by the 
merger are clearly visible. 

4.3. Distorted velocity field 

To heavily disturb the velocity field of a spiral galaxy, a 
tidal interaction with a galaxy of similar or higher mass 
is necessary. In this section we use the velocity field of 
the small galaxy B in simulation 2 after its first passage 
through galaxy A, i.e. about 200 Myr after the first en- 
counter. 

In this case the VF shows no regular pattern of rota- 
tion anymore. Both, at low and intermediate redshift the 
irregularity is visible, although at z = 0.5 the substruc- 
ture is completely smeared out leaving a single 'velocity 
blob' (see Fig. [T5)l . The bad resolution of only 5x5 pix- 
els at this redshift together with the seeing is the reason 
for that. Also the better resolution of 0.25" does not im- 
prove this image significantly. Obviously, such a galaxy 
would not be used for a TuUy-Fisher study and could not 
be mistaken as undisturbed. Also the plots of the radial 
behaviour of the kinemetry parameters clearly reflect the 
distortions in the VF. Due to the low number of pixels 
at intermediate redshifts we constructed this plot only for 
low redshifts. Fig. [16] shows the kinematic parameters at 
z = 0.1. A significant change of the position angle F coin- 
cides with a huge peak in the ratio k^/ki of about 3. Such 
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Fig. 13. The 2D velocity field of a model galaxy A dis- 
turbed by an equal mass merger, 100 Myr after the first 
passage of the two discs is shown in the upper panel. 
The lower panel shows the radial profiles of the kinemet- 
ric properties, calculated using the kinemetry programme. 
The position angle F and the flattening q of the best fit- 
ting ellipses as well as the first and the fifth order Fourier 
terms ki and fcs are plotted as a function of radius. Both 
figures correspond to redshift z=0.5. 



a correlation between the position angle and k^/ki was 
also reported and discussed by Krajnovic et al. (2006). 

Also for this interaction scenario the question arises for 
how long the 2D velocity field maintains the severely dis- 
torted shape. We found that in this case the distortions of 
the velocity field are so strong, that the galaxy at no point 
between the first passage and the final merging could be 
classified as 'disturbed' or even 'regular' anymore. Before 
the first pericentre passage the galaxy features an undis- 
turbed rotation, which is increasingly disturbed when the 
tidal forces of galaxy A become stronger. Once galaxy B 
passes the Milky Way type galaxy the velocity field is 
heavily distorted, as discussed above. 

4.4. Discussion of observational biases 

Throughout the paper we have assumed that the velocity 
field of the galaxy is always observable out to four disc 
scale lengths, which in real observations is of course due 
to limited sensitivity not always possible. Therefore, de- 
pending on the instrument, one might just observe the 
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Fig. 14. 2D velocity field of a modelled disc galaxy dis- 
turbed by an ongoing minor merger event as seen at red- 
shift 2 using SINFONI at VLT. Overlayed as contours are 
the best fitting ellipses from the kinemetry analysis. 




Fig. 16. Radial profiles of the kincmetric properties, cal- 
culated using the kinemetry programme for galaxy B 200 
Myr after a major merger with galaxy A for redshift 
z = 0.1. The position angle F and the flattening q of the 
best fitting ellipses as well as the first and the fifth order 
Fourier terms k\ and fcs are plotted as a function of radius. 




Fig. 15. 2D velocity field of a model galaxy B heavily dis- 
turbed by a merger with galaxy A (200 Myr after the first 
passage through the disc of galaxy A) shown at redshift 
0.1 (left) and 0.5 (right). 



inner parts of the VF. As a consequence distortions in the 
outskirts can be missed. Here we were interested in the 
principle effects of limited resolution on the appearance of 
the 2D velocity fields, without specifying a special instru- 
ment, i.e. for an idealized observation. In a future work 
we want to study flux-limited velocity fields for special 
instruments for comparison with observed data. 

A second parameter that changes the appearance of 
the 2D velocity field is the inclination of the galaxy. 
Throughout the paper we have adopted an inclination 
i=35°. For an undisturbed galaxy typically a correction 
with the sine of the inclination is applied to the rotation 
curve. The appearance of a regular 2D field is not severely 
affected by inclination. In the case of an interacting galaxy, 
however, the appearance of the velocity field changes with 
inclination. A systematic investigation of the effects of in- 
clination is, however, difficult, as they depend on the spe- 



cific interaction geometry. We sketch here the influence of 
inclination on the VF presented in Sect. 14.21 Instead of 
i=35°, we now choose i=80°, i.e. nearly edge-on. The un- 
equal mass edge-on merger of simulation 2 produces VF 
distortions mainly in the plane of the disc. These distor- 
tions are therefore well visible at low inclinations, see Fig. 
|4] but get less prominent for higher inclinations. In Fig. 
[TTI we show the same snapshot as in Fig. |4] but at an in- 
clination of i=80° instead of i=35°. At ^ = 0.5 the VF 
even appears regular, hence the interaction might not be 
recognized. 

4.5. Implications for distant Tully-Fisher studies 

Many recent studies have analysed the luminosity evolu- 
tion of galaxies via the Tully-Fisher relation (e.g. Ziegler et 
al. 2003, Bohm et al. 2004, Bamford et al. 2005, Nakamura 
et al. 2006, Flores et al. 2006). Some of these studies also 
focussed on environmental effects, i.e. differences in the 
Tully-Fisher relation between field and cluster population. 
Ziegler et al. (2003) and Nakamura et al. (2006) find no 
significant differences between field and cluster galaxies, 
whereas Bamford et al. (2005) claim that galaxies in clus- 
ter are on average brighter than their field counterparts. 
The Tully-Fisher relation (TFR) is also extensively used 
as test-bench for galaxy formation models. Portinari & 
Sommer-Larsen (2007) investigated the redshift evolution 
of the Tully-Fisher relation in cosmological simulations. 
They find an offset between the observed and simulated 
Tully-Fisher relation at z=0. The evolution they find is in- 
termediate between diverse observational results. Parts of 
the discrepancies between the various observational results 
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Fig. 17. 2D velocity field of a modelled disc galaxy dis- 
turbed by an ongoing minor merger event at redshift 0.1 
(top) and 0.5 (bottom). In contrast to Fig. [4] the inclina- 
tion in this figure is i=80°. 



may be attributed to the way distortions in the velocity 
fields of the galaxies are accounted for. 

The results from this work suggest that galaxies from 
an observed sample which are misclassified as undisturbed 
can introduce a systematic error on the TuUy-Fisher di- 
agram. A TuUy-Fisher analysis assumes that the galaxies 
are in virial equilibrium, which is revealed by a smooth, 
regular and symmetric rotation curve, that rises rapidly 
in the inner part and turns over into a flat part. This 
"flat" or maximum circular velocity is one parameter of 
the TFR, the other being the luminosity. In a disturbed 
system, the measured Ymax is not an accurate estimate for 
the virial velocity of the halo anymore and, hence, no good 
proxy for the total mass of the system. Also, the amount 
of luminosity of the galaxy is influenced by the interaction 
(see Kapferer et al. 2005). To which extent both quantities 
are affected by the interaction depends on the interaction 
geometry and the time. 

Therefore, if distortions or peculiarities are visible in 
the RC, the galaxy would not be included in a TF analysis. 
However in this study, we have shown that at intermediate 
redshifts due to lack of resolution, signatures of an inter- 
action are smeared out, so that a distorted galaxy mimics 
a regular one in several cases. If such disturbed (but un- 
recognized) systems are used or a TuUy-Fisher study, the 
scatter will be significantly increased. In contrast to slit 
observations of a one-dimensional RC, the availability of 
a two-dimensional VF allows a much better assessment 
of irregularity. Nevertheless, in case of a small and faint 



galaxy, where the VF spreads over a few pixels only, an 
interaction could be missed as well. The inclusion of such 
intrinsically disturbed objects could therefore mimic an 
evolution at the faint end of the TuUy-Fisher relation. 

We show in Fig. [TH] how a galaxy moves in the TuUy- 
Fisher diagram if Ymax values are used that do not cor- 
respond to the virial velocity. First, we place the undis- 
turbed galaxies onto the TFR under the assumption of 
a constant mass-to-light ratio. The simulated Milky Way 
galaxy A and the less massive galaxy B then follow a lin- 
ear relation with slope of 3, i.e. L oc Vf^^^. This is in 
good agreement with observed local values. Pierce & TuUy 
(1992) get a value of 2.99 for their sample, other authors 
get values between 2.96 and 3.2 (see Ziegler et al. 2002 and 
references therein). In Fig.[T8]we then sketch how interac- 
tions change the positions of the galaxies. Note that these 
trends are taken from the interaction scenarios studied 
in this work. Other interaction geometries or times might 
lead to a different effect. The Milky Way type galaxy A 
moves to the right, if the minor merger presented in Sect. 
14.21 is not recognized. As the interaction also triggers star 
formation, the luminosity increases at the same time. In 
the case studied here, the effect on Ymax is greater than 
the brightening of the galaxy. For the small galaxy B in the 
case of a major merger the interactions lead to an under- 
estimation of the true Ymax- At the same time, the lumi- 
nosity stays roughly constant or even decreases by a small 
amount, as stars are lost by the merger. Assuming this 
trend holds more generally and that distortions are more 
easily missed at higher redshift for small galaxies, this ef- 
fect could mimic a luminosity evolution at the faint end 
of the TuUy-Fisher relation. Therefore, a careful analysis 
is necessary to prevent the inclusion of disturbed systems 
in distant TuUy-Fisher studies. 

5. Summary and conclusions 

We have investigated 2D velocity fields of isolated and in- 
teracting spiral galaxies using N-body/SPH simulations. 
We focussed on the question how the full 2D velocity field 
of a galaxy can be used to gain information on its inter- 
nal kinematics. This issue was analysed with special em- 
phasis on distant TuUy-Fisher studies. To summarise and 
conclude: 

— We found that with the help of 2D velocity fields the 
nature of the interaction becomes a lot more accessible 
than with rotation curves from simple long-slit spec- 
troscopy. Quantitative analysis such as an harmonic 
expansion, which is, for example, used by the kineme- 
try package of Krajnovic et al. (2006), offer an addi- 
tional possibility to identify distortions in a VF. 

— Tidal interactions lead to an asymmetric velocity field, 
where the side pointing at the interaction is disturbed 
while the side remote from the interaction stays rel- 
atively unaffected. This behaviour can also be found 
in the rotation curve of the system (Kronberger et al. 
2006). 
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Fig. 18. Position of the two undisturbed model galaxies in 
the luminosity-rotational velocity plane (red cross: galaxy 
A; blue asterisk: galaxy B). The two galaxies lie on a TuUy- 
Fisher relation with slope 3 (dotted line), consistent with 
observations of local galaxy samples. The arrows indicate 
the trends in which way the interactions studied in this 
work change the positions of the galaxies. For the Milky 
Way type galaxy, the investigated interaction is the minor 
merger with galaxy B, while for model galaxy B it is the 
major merger case with a second galaxy B. These trends 
can be observed at all redshifts but in a different strength. 



misclassification as 'undisturbed rotation' is possible 
at higher redshifts. 
— With the assumption of a constant mass-to-light ra- 
tio, we showed that the undisturbed model galaxies lie 
on a TuUy-Fisher relation with a slope of 3, which 
is consistent with observations of local galaxy sam- 
ples. Using systems with unrecognized distortions for a 
Tully-Fisher study will significantly increase the scat- 
ter in the observed relation. As it is more likely, that 
an interaction is missed for small, faint galaxies, the 
inclusion of disturbed systems can mimic an evolution 
in the faint end of the Tully-Fisher relation. 
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— By shifting the VF artificially to higher redshifts we 
found that, although small scale structures in the VF 
are smeared out, distortions are still visible at interme- 
diate redshifts for large galaxies. In the case of small 
galaxies even strong distortions are not visible in the 
velocity field at z « 0.5 with currently available angu- 
lar resolution. 

— Severely distorted kinematics, i.e. velocity fields with 
no or just a small rotational component are identifiable 
also for small galaxies at intermediate redshift. Thus 
they should not be mistaken as undisturbed. 

— If the galaxy is undisturbed the quantities derived from 
the velocity field, as e.g. the maximum rotational ve- 
locity Vmaa;, do not show Systematic variations with 
resolution and are therefore useable for Tully-Fisher 
studies. 

— We showed that with the help of adaptive optics near- 
infrared spectrographs (e.g. SINFONI) it is possible to 
study the internal kinematics of galaxies even at high- 
redshift (z^2), given that the observed flux is sufficient 
to construct a 2D velocity field. Ongoing merger events 
are then clearly visible in the VF. 

— Disturbed velocity fields settle again to a relaxed state 
after about 1 Gyr. However, already after several hun- 
dred Myr the distortions can get damped such that a 
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